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Throughout physics, stable composite objects are usually formed via attractive forces, which allow the con- 
stituents to lower their energy by binding together. Repulsive forces separate particles in free space. However, 
in a structured environment such as a periodic potential and in the absence of dissipation, stable composite ob- 
jects can exist even for repulsive interactions. Here we report on the first observation of such an exotic bound 
state, comprised of a pair of ultracold atoms in an optical lattice. Consistent with our theoretical analysis, these 
repulsively bound pairs exhibit long lifetimes, even under collisions with one another. Signatures of the pairs are 
also recognised in the characteristic momentum distribution and through spectroscopic measurements. There 
is no analogue in traditional condensed matter systems of such repulsively bound pairs, due to the presence of 
strong decay channels. These results exemplify on a new level the strong correspondence between the optical 
lattice physics of ultracold bosonic atoms and the Bose-Hubbard model0,01, a correspondence which is vital 
for future applications of these systems to the study of strongly correlated condensed matter systems and to 
quantum information. 
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Cold atoms loaded into a 3D optical lattice provide a reali- 
sation of a quantum lattice gas| 1, 2]. An optical lattice can be 
generated by pairs of counterpropagating laser beams, where 
the resulting standing wave intensity pattern forms a periodic 
array of microtraps for the cold atoms, with period a given by 
half the wavelength of the light, A /2. The periodicity of the 
potential gives rise to a bandstructure for the atom dynamics 
with Bloch bands separated by band gaps, which can be con- 
trolled via the laser parameters and beam configuration. The 
dynamics of ultracold atoms loaded into the lowest band of a 
sufficiently deep optical lattice is well described by the Bose- 
Hubbard model with Hamiltonian| l||3J 



6 = -J £ b]bj + 1 £ b\bi (b]h mBJh 



(i) 



Here bi (bj) are destruction (creation) operators for the 
bosonic atoms at site i, J/h denotes the nearest neighbour 
tunnelling rate, U the on-site collisional energy shift, and 
e, the background potential. The high degree of control 
available over the parameters in this system, e.g., changing 
the relative values of U and J by varying the lattice depth, 
Vq, has led to seminal experiments on strongly correlated 
gases in optical lattices, e.g., the study of the superfluid- 
Mott insulator transitionfl, the realisation of ID quantum 
liquids with atomic gases^ El ( see also0, |H), and the 
investigation of disordered systems|9]. 3D optical lattices 
have also opened new avenues in cold collision physics and 

chemistry E3 E3 Q • 

A striking prediction of the Bose-Hubbard Hamiltonian Q 
is the existence of stable repulsively bound atom pairs. These 
are most intuitively understood for strong repulsive interaction 
| U | 3> J, U > 0, where an example of such a pair is a state 

of two atoms occupying a single site, |2,) = (»T 2 |vac))/-\/2- 
This state has a potential energy offset U with respect to states 



where the atoms are separated (see Fig. la). The pair is un- 
able to decay by converting the potential energy into kinetic 
energy, as the Bloch band allows a maximum kinetic energy 
for two atoms given by 8/, twice its width. The pair can 
move around the lattice, with both atoms tunnelling to a neigh- 
bouring site (Fig. lb), but the atoms cannot move indepen- 
dently. The stability of repulsively bound pairs is intimately 
connected with the absence of dissipation, e.g., in contrast to 
solid state lattices, where interactions with phonons typically 
lead to rapid relaxation. 

We obtain experimental evidence for repulsively bound 
pairs with a sample of ultracold 87 Rb atoms in a 3D optical 
lattice with lattice period a = 415.22 nm. The key tool used 
to prepare and observe the pairs is their adiabatic conversion 
into chemically bound dimers using a magnetic-field sweep 
across a Feshbach resonance 1141 Il5l lie* Il7l Il8l Il9i 12011 (at 
1007.40 G, see 1131 12 1IH . The initial state is prepared from 
a pure sample of Rb2 Feshbach molecules in the vibrational 
ground state of the lattice where each lattice site is occupied 
by not more than a single molecule| 13]. Sweeping across 
the Feshbach resonance we adiabatically dissociate the dimers 
and obtain a lattice correspondingly filled with 2 x 10 4 atom 
pairs, at an effective filling factor of typically 0.3. In all our 
experiments reported here, the preparation takes place at a 
deep lattice of Vo = 35 E r (E r — 2n 2 h 2 /mX 2 , where m is the 
mass of the atoms), corresponding to J/h «2ix 0.7 Hz and 
U/J = 3700. Away from the Feshbach resonance, the effec- 
tive interaction between the atoms is repulsive with scattering 
length a s — +100ao (where «o is me Bohr radius). 

In order to demonstrate the stability of repulsively bound 
pairs, we lower the lattice potential in one direction in 1 ms 
to a depth of Vo = \QE r . This increases dramatically the tun- 
nelling rates along this direction to J /h rj 2n x 63 Hz, poten- 
tially allowing the atom pairs to quickly separate. After a vari- 
able hold time we determine the number of remaining pairs. 
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FIG. 1 : Atom pairs in an optical lattice, (a) Repulsive interaction 
(scattering length a > 0) between two atoms sharing a lattice site in 
the lowest band (n = 0) gives rise to an interaction energy U. Break- 
ing up of the pair is suppressed due to the lattice band structure and 
energy conservation, (b) The pair is a composite object that can tun- 
nel through the lattice, (c) Long lifetime of repulsively bound atom 
pairs which are held in a 3D optical lattice. The potential depth is 
(10±0.5)£ r in one direction and (35 ± 1.5) E r in the perpendicular 
directions. Shown is the remaining fraction of pairs for a scattering 
length of 100 ao (open diamonds) and a scattering length of about 
(0 ± 10) ao (filled circles) as a function of the hold time. The lines 
are fit curves of an exponential (dashed line) and the sum of two ex- 
ponentials (solid line). 



This is done by adiabatically raising the lattice to its full ini- 
tial depth of Vq = 35 E r , and converting doubly occupied sites 
to Feshbach molecules with near unit efficiency |13J. A purifi- 
cation pulse 1 13] then removes all remaining atoms due to dis- 
sociated pairs. Afterwards the molecules are again converted 
back into atoms, and can then be detected via conventional 
absorption imaging. 

The results of these lifetime measurements are shown in 
Fig. lc. For repulsive interaction (a s = 100ao) the atom 
pair sample exhibits the remarkably long lifetime of 700 ms 
(exponential fit). This lifetime is mainly limited by inelastic 
scattering of lattice photons llOI . and greatly exceeds the cal- 
culated time for an atom to tunnel from one site to the next, 
2nh/(4J) ~ 4ms. In contrast, if we turn off the on-site inter- 
action by tuning the scattering length near zero, we observe 
a much faster decay in the number of doubly occupied sites 
due to the rapid diffusion of unbound atoms through the lat- 
tice (see Fig. lc). This observation clearly demonstrates that 
the stability of the pairs is induced by the on-site interaction 
U. 

We can more deeply understand these repulsively bound 
pairs through an exact solution of the two particle Lippmann- 
Schwinger scattering equation based on the Bose-Hubbard 
model. We write the two atom wavefunction as *P(x,y), 
where the positions of the two particles are denoted x = 
Y,i x &i an d y = T^iyfih with e i being the primitive lattice vec- 
tors, and xi,yi integer numbers. Introducing centre of mass, 
R = (x + y)/2, and relative coordinates, r = x — y, we can 
solve the Schrodinger equation with the ansatz ^(Xjjr) = 
exp(/KR)i//)r(r), where K is the quasi-momentum of the cen- 
tre of mass motion. We derive two types of solutions (for 
details see appendix lAl. each of which are eigenstates of the 
centre of mass quasi-momentum K. These states, as illus- 
trated in Fig. 2a, correspond to (i) unbound scattering solu- 



tions (shaded area in Fig. 2a), where the two particles move 
on the lattice, and scatter from each other according to the in- 
teraction U, and (ii) repulsively bound pairs for which the pair 
wavefunction Vk"(r) is square integrable. In ID and 2D, states 
of repulsively bound pairs always exist for non-zero U, while 
in 3D they exist above a critical value t/ cr it ~ 0.57. 

In this letter, we primarily focus on the ID situation, which 
in the experiment corresponds to a low depth of the lattice 
along one direction, whilst the lattice in the perpendicular 
directions remains very deep (35E r ). Here the energy of 

the bound pairs is E(K) 
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This is plotted in Fig. 2a as the Bloch band of a stable com- 
posite object above the continuum of two particle scattering 
states. In the limit of strong interaction, U » 7, this re- 
duces to E(K) ~ 47 + U + (AJ 2 /U ) (1 + cosKa) , which shows 
that the bound pairs indeed have binding energy w U and hop 
through the lattice with an effective tunnelling rate J 2 / (hU). 

Figure 2b shows the pair wavefunctions y/jr(r) for repul- 
sively bound pairs (a s = 100ao) in ID with K = 0, for U/J = 
30 (Vb = 10£ r ) and U/J = 3 (V = 3E r ). For large U/J, 
bound pairs essentially consist of two atoms occupying the 
same site, whereas for small U /J, the pair is delocalised over 
several lattice sites. The corresponding quasi-momentum dis- 
tribution can be found from the Fourier transform \j/a(k) of 
the pair wavefunction (see Fig. 2c). Since K = 0, |y/o(&)| 2 
is also equal to the single particle quasi-momentum distribu- 
tion. When the two particles are localised on the same site, the 
quasi-momentum distribution is essentially flat. However, fol- 
lower U/J the wavefunction is characteristically peaked at the 
edges of the Brillouin zone. This occurs because the energy 
of the repulsively bound state is above that of the continuum, 
and thus the contribution to the corresponding wavepacket of 
single particle quasi-momentum states with higher energy is 
favoured. In contrast, if we had U < 0, the pair would be 
attractively bound, and would have energy lower than that 
in the continuum. Thus contributions from the low energy 
quasi-momentum states would be favoured, leading to a sin- 
gle peak in the centre of the Brillouin zone. In both cases, 
the amplitude of the peaks grows with increasing width 47 of 
the Bloch band. In general, the stable bound pairs will not 
be prepared in a fixed quasi-momentum state K in an exper- 
iment, but rather in a superposition of different momentum 
states. For non-zero K, the peaks in the single particle quasi- 
momentum distribution are translated by K, but their strength 
is also reduced. As a consequence, for typical symmetric dis- 
tributions of K, the peak at the edge of the Brillouin zone re- 
mains present, however, is less strong than the optimal case of 
vanishing K. We have verified this using many-body numer- 
ical simulations, which are performed using time dep endent 
density-matrix renormalisation group methods[22, 23l l24ll . 

We have experimentally investigated the quasi-momentum 
distribution of the pairs in various regimes by mapping it onto 
a spatial distribution, which we measured using standard ab- 
sorption imaging. For this, we first adiabatically lower the 
lattice depth in the X - direction (see Fig. 3a) to a pre- 
chosen height while the lattice depth in the other two direc- 
tions are kept high (35 E r ). This will prepare repulsively 
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FIG. 2: Atom pair states in ID. (a) Energy spectrum of the ID Hamiltonian for U/J = 5 (Vo ~ 4E r ) as a function of centre of mass quasi- 
momentum K. The Bloch band for repulsively bound pairs is located above the continuum of unbound states. The grey level for the shading 
of the continuum is proportional to the density of states, (b) The pair wavefunction Vb( r )> showing the amplitude at each site with U/J — 30 
(Vo f» lOE r , narrow bars) and U/J = 3 (Vo ~ 3E r , wide bars), (c) The square modulus of the corresponding momentum space wavefunctions 
|y/o(fc)| 2 , which are equivalent to the single particle momentum distributions, since K = 0. Note the characteristic peaks at the edge of the 
Brillouin zone. 



bound pairs at the chosen lattice depth. We then turn off the 
lattice rapidly enough so that the pair wavefunction cannot 
change, but slowly with respect to the bandgap, so that single- 
particle quasi-momenta are mapped to real momenta 1 25l 12611 . 
We have typically employed linear ramps with rates of 0.2 
E r /fis. The resulting momentum distribution is converted to 
a spatial distribution after ~ 15 ms time of flight. 

Figure 3 (a-c) shows typical measured quasi-momentum 
distributions which were obtained after adiabatically lowering 
the lattice depth in the X - direction to lowest lattice depths 
below 3E r . If only empty sites and sites with single atoms 
are present in the lattice, then the first Brillouin zone is ho- 
mogenously filled [25] (Fig. 3 a). For repulsively bound pairs 
the momentum distribution is, in general, peaked at the edges 
of the first Brillouin zone (Fig. 3 b), whereas for attractively 
bound pairs it is peaked in the center of the first Brillouin zone 
(Fig. 3 c). In order to change the interaction between the 
atoms from repulsive to attractive, we change the scattering 
length making use of the Feshbach resonance|21] at 1007.40 
G. Figures |3]d and e show the dependence on lattice depth 
Vo of the single particle quasi-momentum distribution for re- 
pulsively bound pairs from experiment and numerical simu- 
lation, respectively. As expected, the peak structure is more 
pronounced for lower values of Vo, and diminishes for larger 
Vo. This characteristic is a clear signature of the pair wave- 
function for repulsively bound pairs. 

We also performed spectroscopic measurements, determin- 
ing the binding energy from dissociation of the pairs by mod- 
ulating the depth of the lattice at a chosen frequency. On res- 
onance, the modulation allows pairs to release their binding 
energy. Figure 4 a shows the number of remaining pairs as a 
function of the modulation frequency. This was repeated for a 
variety of lattice depths Vo in one direction while keeping the 
lattice in the other two directions at 35 E r . The behaviour of 
the binding energy as a function of the lattice depth provides 
an additional key signature of repulsively bound pairs. As 
shown in Fig. 4b, the resonance positions are in good agree- 
ment with numerical simulations and essentially coincide with 
interaction energy, U . 

It is important to note that for sufficiently large U/J, repul- 
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FIG. 3: Quasi-momentum distribution of atoms in the lattice. Im- 
ages (a) - (c) show absorption images of the atomic distribution after 
release from the 3D lattice and a subsequent 15 ms time of flight. The 
horizontal and vertical black lines enclose the first Brillouin zone, (a) 
Distribution when lattice sites are occupied by single atoms, (b) Dis- 
tribution for repulsively bound atom pairs (see text for details), (c) 
Same as (b) but pairs are attractively bound. Panels (d) [experiment] 
and (e) [numerical calculation] depict the quasi-momentum distribu- 
tion for pairs in the X-direction as a function of lattice depth Vo, a f ter 
integration over the F-direction. 



sively bound pairs are stable under collisions with each other. 
This is particularly evident in the limit U » J where, by 
energy arguments, the elastic scattering between pairs is the 
only open channel. This means that even a relatively dense 
quantum lattice gas of these objects can be long-lived. When 
the lattice height is lowered so that U /J becomes sufficiently 
small, it is possible for a certain fraction of the pairs to dis- 
sociate by collision with other pairs. In our experiments, we 
observe the onset of this behaviour for lattice depths lower 
than 6E r , i.e., U/J ~ 9. The dynamics of the collisions and 
details of the decay depend crucially on lattice depth and the 
local density of pairs across the lattice. Further details of these 
processes will be discussed elsewhere. 

In conclusion, we have demonstrated the formation of a 
novel composite object in an optical lattice: a stable bound 
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FIG. 4: Modulation spectroscopy of repulsively bound pairs, (a) Typ- 
ical resonance dip showing the remaining number of atom pairs as a 
function of the modulation frequency, for Vn ~ 6E r . The black line 
is a Gaussian fit, a choice which was justified by numerical calcula- 
tions, (b) Plot showing the measured resonance frequencies (filled 
circles) as a function of the lattice depth. They show good agreement 
with numerical simulations (crosses) and also coincide with the on- 
site collisional energy shift U (line). Experimental error bars corre- 
spond to the 95% confidence interval for the Gaussian fit parameters 
of the resonance dips. 



state that arises from the lattice band structure and repulsion 
between the constituents. Although no direct analog to re- 
pulsively bound atomic pairs is known to exist, the formation 
of a metastable state is reminiscent of trapping light in pho- 
tonic bandgap materials 112711 . or extended lifetimes of excited 
atoms in cavity quantum electrodynamics 1 28]. In both cases, 
decay is suppressed by restriction of the accessible light field 
modes. Stable repulsively bound objects should be viewed 
as a general phenomenon and their existence will be ubiq- 
uitous in cold atoms lattice physics. They give rise also to 
new potential composites with fermions|29] or Bose-Fermi 
mixturesl 30], and can be formed in an analogous manner with 
more than two particles. The stability of repulsively bound ob- 
jects could thus be the basis of a wealth of new quantum many 
body states or phases. 



tor AK¥(x)=Ef =1 cos(Ke;/2) ^(x+eO+^x-e/) -2¥(x)] 
denotes the discrete lattice Laplacian with d the dimensional- 
ity in the cubic lattice, and 5x. y is a Kronecker delta. Writ- 
ing the wavefunction in relative and centre of mass coordi- 
nates *P(x,y) = exp(/KR)y/R - (r), the Schrodinger equation 
(IA 1 1 then reduces to a single particle problem in the relative 
coordinate 



[-2JA? + E K + U8 r .o] Wk{v) = £ft(r) (A2) 
with E K = 4JY,i [1— cos(Ke,/2)] being the kinetic energy of 
the centre of mass motion. 

The short range character of the interaction potential allows 
for a resummation of the perturbation expansion generated by 
the corresponding Lippmann-Schwinger equation. We obtain 
the scattering states 



V/ (+) (r) =exp(iTff)-8Jr// £ (K)G K (£,r) (A3) 



with scattering amplitude 



/*(K) 



1 



U/(2J) 



\% l-G K (E,0)U 



(A4) 



where the total energy E = e kK + E K , and £ kjK = 
4/£ I cos(Ke,/2) [1 — cos(ke,-)]. Furthermore, G K {E,r) de- 
notes the Greens function of the non-interacting problem, 
which in Fourier space takes the form Gn(E,k) = \/{E — 
£k, K + n7)- 

The energy spectrum for these states in one dimension is 
shown as a function of K by the shaded region in Fig. 2a. 
In addition, the pole in the scattering amplitude indicates the 
presence of an additional bound state. The energy E bs of the 
bound state is determined by G K (E bs ,Q)U = 1 and the bound 
state wavefunction takes the form \j/ hs (r) = cG%(E bs ,r) with c 
being a normalisation factor. 



APPENDIX A: EXACT SOLUTION FOR SINGLE PAIR 
BOUND STATE 

Within the Bose-Hubbard Model (1) the Schrodinger equa- 
tion describing two particles in a homogenous optical lattice 
takes the form 

[-J (A° + A°) + U8^\ *F(x,y) = E *P(x, y), (Al ) 

where the vectors x and y describe the positions of the 
two particles as defined in the main text. The opera- 
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